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^ . Abstract. 

We review some recent developments in theoretical studies on the connection between the 
|_ progenitor systems of Type la supernovae (SNe la) and the explosion mechanisms. (1) DD- 
subCh: In the merging of double C+O white dwarfs (DD scenario), if the carbon detonation 
is induced near the white dwarf (WD) surface in the early dynamical phase, it could result in 
the (effectively) sub-Chandrasekhar mass explosion. (2) DD-Ch: If no surface C-detonation is 
ignited, the WD could grow until the Chandrasekhar mass is reached, but the outcome depends 
on whether the quiescent carbon shell burning is ignited and burns C+O into O+Ne+Mg. (3) 
SD-subCh: In the single degenerate (SD) scenario, if the He shell-flashes grow strong to induce 
a He detonation, it leads to the sub-Chandra explosion. (4) SD-Ch: If the He-shell flashes are 
not strong enough, they still produce interesting amount of Si and S near the surface of C+O 
WD before the explosion. In the Chandra mass explosion, the central density is high enough 
to produce electron capture elements, e.g., stable 5S Ni. Observations of the emission lines of Ni 
( in the nebular spectra provides useful diagnostics of the sub-Chandra vs. Chandra issue. The 
recent observations of relatively low velocity carbon near the surface of SNe la provide also 
interesting constraint on the explosion models. 
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1. Introduction 

The observed features of Type la supernovae (SNe la) have been well-understood by a 
thermonuclear explosion of a carbon-oxygen (C+O) white dwarf (WD). Both the Chan- 
drasekhar mass [Chandra (Ch) model] and the sub-Chandrasekhar mass [sub-Chandra 
(subCh) model] have been presented (e.g., ILivio 20001 - However, there has been no 
clear observational indication as to how the WD mass grows until carbon ignition; i.e., 
whether the WD accretes H/He-rich matter from its binary companion [single degen- 
erate (SD) scenario], or two C+O WDs merge [douWe degenerate (DD) scenario] (e.g., 
|Hillebrandt fc Niemeyer 2000[ INomoto et d. 19971 2000, 2009, lArnett 1996) . Even be- 
fore these issues are resolved, several candidates of super- Chan drasekhar mass explosions 
have been observed (e.g.. IHachisu et al. 2012l and references therein). 

Recent modeling shows that DD merging could result in both the Chandra and the 
(effectively) sub-Chandra explosions. The SD scenario could also result in both the Chan- 
dra and sub-Chandra explosions. Here we review such cross connections between (DD, 
SD) scenarios and (Chandra, sub-Chandra) models, and some observational constraints. 

(1) DD-subCh: In the DD scenario, if the carbon detonation is induced near the WD 
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surface in the early dynamical phase, it could result in the (effectively) sub-Chandra 
explosion. (2) DD-Ch: If no detonation is induced, the WD could grow until the Chan- 
drasekhar mass is reached. The outcome depends on whether the quiescent carbon shell- 
burning is ignited and burns interior C+O into O+Ne+Mg. (3) SD-subCh: In the SD 
scenario, if the He shell-flashes grow strong to induce a He detonation, it leads to the 
sub-Chandra explosion. (4) SD-Ch: If the He-shell flashes are not strong enough to induce 
a He detonation, such flashes produce interesting amount of intermediate mass elements, 
including Si and S, as unburned material near the surface of C+O WD. 

2. C+O Double Degenerates to Sub-Chandra Mass Explosion 

In the DD scenario fllben fc Tutukov 19841 IWebbink 1984|) . two C+O WDs form a 
close binary system after the common envelope phase and get closer by losing orbital 
angular momentum due to gravitational radiation. Eventually, the less massive WD with 
a mass M2 fills its Roche lobe and dynamically evolves into the formation of a massive 
disk/envelope around the more massive WD with a mass M\. 

Such a dynamical evolution of WD binary has been studied by number of authors 
(jBenz et al. 19901|Segretain et al. 19971 IGuerrero et al. 2004l|Shioya et al. 2007[lPakmor et al. 2010]) 
using the Smoothed Particle Hydrodynamics (SPH) method. 

Surface Carbon Detonation: The important question is whether the merging pro- 
cess ignites a surface carbon detonation in the WD. During the early merging process, the 
shock heating increases the temperature in the surface C+O layer of the primary (i.e., 
more massive) WD. If the temperature becomes high enough to induce C-detonation 
(jPakmor et al. 20101 adopted the critical temperature of 2.5 xlO 9 K), the detonation 
wave propagates through the central region of the primary WD. Eventually, the whole 
primary WD is detonated and disrupted. 

This model is essentially the sub-Chandra mass explosion, because the detonated WD 
has a sub-Chandra mass of Mi = O.9-1.1M and the central density is as low as ~ 10 7 
g cm~ 3 (Pa kmor et al. 20101 2011). The explosion produce a larger amount of 56 Ni for 
larger M\. 

3. Double Degenerates to Chandra Mass Explosion or Collapse 

Suppose the surface C-detonation is not triggered in the early dynamical phase of merg- 
ing. Then the next important question is whether the merging process ignites not the C- 
detonation but "quiescent" off-center carbon burning in the WD. Once carbon is ignited, 
carbon flame moves inward by heat conduction to reach the center (jSaio fc Nomoto 1985|, 
1998). The released nuclear energy is lost in neutrinos, and the C+O WD is converted 
into the O+Ne+Mg WD non-explosively. The ONeMg WDs eventually collapse due to 
electron capture rather than exploding as SNe la (jNomoto 1984) 1987). 

3.1. Carbon Shell Burning and Chandra Explosion or Collapse 

Here the carbon ignition temperature T- lgn is defined by ec+c = e i/, where ec+c and e„ 
denote the nuclear energy generation rate and the neutrino energy loss rate, respectively. 
At p ~ 1-3 x 10 6 g cm" 3 , T ign ~ 6 x 10 8 K (jNomoto fc Iben 19851 IKawai et a/.TJ87l 
lYoon et al.~2007\i . 

For T > Ti gn , ec+c > £v and the carbon flash is ignited, and the conductive carbon 
flame propagates inward. If T < T| gn , on the contrary, ec+c < and the neutrino 
cooling dominates to induce the gradual contraction of the C-rich envelope. Then carbon 
flame is not formed. 
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Figure 1. The local peak temperature T pea k attained 2.5 min after the start of merging of 
double C+O white dwarfs with masses (Mi, M 2 ) (|Nakasato fc Nomoto 2011[) . For T pcak > T ign , 
off-center carbon burning is ignited. 



There are two possible cases for the off-center carbon ignition to occur. 
Case 1: When the falling materials from the less massive star compress the outer layer 
of the more massive star, the material could be heated up to high enough temperature 
to ignite carbon burning. 

Case 2: Later, if the accretion rate exceeds a critical rate of 2.7 x 1O _6 M0 yr _1 , compres- 
sional heating exceeds radiative cooling and leads to carbon ignition (N omoto fc Iben"l 985, 
IKawai et al. 1987|) . 

For Case 1, Yoon et al. (2007) calculated the merging process until 5 min after its start 
and showed that a quasi-static equilibrium configuration is reached consisting of a cold 
core, hot envelope, and a disk. The peak temperature T poa k reaches a stationary value. 
For Mi = 0.9M© and M 2 = O.6M , T pca k is marginally lower than T ign . Whether the 
later off-center carbon ignition of Case 2 takes place depends on the effective accretion 
rate, which needs further study (e.g., |Shioya et al. 2007|IShen et al. 2012|) . 

For Case 1, Nakasato & Nomoto (2011) have recently conducted the SPH simula- 
tion (|Nakasato fc Nomoto 2003j) with the number of particles N = 300,000 and N = 
1,000,000 for various combinations of (Mi, M 2 ). The artificial viscosity is treated to 
minimize numerical effects according to the hybrid scheme proposed by Rosswog et al. 
(2010). The Helmholtz equation of state is used ( |Timmes fc Swesty 2000[ ), and the WD 
is assumed to be composed of 50 % of carbon and 50 % of oxygen. 

The local peak temperatures shown in Figure Q] are obtained from particles at p ~ 1- 
3 x 10 6 g cm~ 3 at 2.5 min after the start of merging. T pea k has already reached its 
stationary value, which is confirmed with some test runs calculated until 10 min after 
the merging. It is seen that T pea k is determined mainly by Mi with small dependence on 
M2- This implies that the gravitational potential of the more massive WD is the main 
factor to determine T pea k. 

Figure [1] shows that T pea k > Ti gn for most cases of (Mi, M 2 ). Thus carbon flash 
will take place to form a carbon flame that propagates inward to convert C+O into 
O+Ne+Mg. 

How the above results depend on the set-up of the merging calculations (jDan et al. 2011|) ? 
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Figure 2. Abundance distribution of the products of He shell flashes in an accreting WD just 
before the explosion (Kamiya & Nomoto 2011). 

The initial separation is set to be a = (0.9i?2)/ r L where R2 is the radius of the less mas- 
sive WD and 7% is the effective Roche lobe radius. Some test runs adopting a larger 
separation confirm that Tpeak is not sensitive to the initial separation. 

Comparisons between the runs with N = 300, 000 and with N = 1, 000, 000 show that 
Tpcak of the higher resolution run is always slightly higher than the low resolution run. 
Thus carbon ignition of Case 1 is quite likely for most cases of (Mx,M2). 

3.2. Further Evolution of Rotating White Dwarfs 

The WD formed from merging must be rapidly rotating. The "SN la mass" of the 
WD with which SN la is triggered is Mr a = 1.48M Q for a uniformly rotating WD 
at the critical rotation (jUenishi et al. 2003)) . This is larger than 1.38M for the non- 
rotating WD. For non-uniform, differentially rotating WDs, Mj a is as large as > 2M 
QYoon fc Langer 2004| IHachisu et al. 2012|) . Therefore, even if no carbon flame is formed, 
the WD may not reach Mj a because of some mass loss after merging. 



4. Single Degenerate to Sub-Chandra Explosion 

In the SD scenario, H-burning produces a thin He layer, and He-flashes are ignited 
when the He mass reaches a certain critical value. The strength depends on the He 
envelope mass M cnv , thus depending on the accretion rate. 

The He envelope mass M env is larger for the slower mass-accumulation rate of the He 
layer Mn e . For A/h > lx 10~ 8 M Q yr _1 , M cnv exceeds a critical value where the density 
at its bottom becomes high enough to induce a He detonation. This would result in the 
sub-Chandra explosion (e.g., S. Sim in this conference). 



5. Single Degenerate to Chandra Mass Explosion 

For higher Mh 6 , the He-shell flashes are not strong enough to induce a He detonation. 
Then, such flashes repeat many times with the increasing WD mass -Mwd toward the 
Chandrasekhar mass. Kamiya & Nomoto (2011) calculated nucleosynthesis in such He 
shell flashes for various set of (Mwd, M env ) (see also lShen fc Bildsten 2007[) . 
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Figure 3. Analysis of the [Ni II] A7378 line profiles in 12 SNe la (|Maeda et al. 2010a|l . The 
velocity is set assuming that the rest wavelength is at 7378 A. (a) Observed line profiles. The 
rest wavelengths of [Fe II] A7155 and [Ni II] A7378 are shown by dotted lines, (b) [Ni II] A7378 
in observations, after removing the underlying continuum (or possible other lines), (c) Synthetic 
line profiles of the [Ni II], depending on the viewing orientation. 

For Miie ~ (0.5 — 1) x 1CT 7 M© yr _1 , the WD is expected to undergo He shell-flashes 
at (M wd/Mq, log(M env /M©)) ~ (1.2,-2.5), (1.3,-3), (1.35,-3.5), (1.38,-4) as M W d 
grows IjKato et al. 2008|> . 

In the early stages of the He shell-flash, the envelope is electron-degenerate and ge- 
ometrically almost flat. Thus the temperature at the bottom of the He-burning shell 
increases because of the almost constant pressure there. Heated by nuclear burning, the 
helium envelope gradually expands, which decreases the pressure. Then, the tempera- 
ture attains its maximum and starts decreasing. The maximum temperature is higher 
for more massive WD and more massive envelope because of higher pressure. 

For higher maximum temperatures, heavier elements, such as 28 Si and 32 S, are syn- 
thesized. However, the maximum temperature is not high enough to produce 40 Ca. After 
the peak, some amount of He remains unburned in the flash and burns into C+O during 
the stable He shell burning. 

In this way, it is possible that interesting amount of intermediate mass elements, in- 
cluding Si and S, already exist in the unburned C+O layer at M r ^ 1.2M©. Such a 
distribution is shown in Figure[5]for Mu a ~ 5 x 10~ 8 M© yr _1 . A part of them might be 
ejected out. Therefore the above quantity of the synthesized elements are overestimated. 

6. Electron Capture in Chandra Mass Models 

Both Chandra and sub-Chandra explosion models can synthesize relevant amount of 
56 Ni for SNe la. However, the amount of other Fe-peak elements differs, because the 
ignition density is different, being as high as > 10 9 g cm -3 in the Chandra model, while 
as low as ~ 10 7 g cm~ 3 in the sub-Chandra model. 

In the Chandra model, the thermonuclear runaway starts with the ignition of defla- 
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Figure 4. Spectrum of SN 2005el at —7 days (solid line), and a matching SYNOW calculation 
(dashed line). (Insets) The spectrum near the location of C II lines (solid line), the SYNOW 
spectra containing only C II (dashed line), and only H I (dotted line). The expected locations of 
C II lines are marked (Folatelli et al. 2012) 

gration (e.g., INomoto et al. 19761 1984). In the high temperature and density bubble, 
materials are incinerated into NSE (nuclear statistical equilibrium) and undergo electron 
capture. Electron capture on free protons and Fe-peak elements leads to the synthesis of 
58 Ni, 54 Fe, and 56 Fe (not via 56 Ni decay). These neutron- rich Fe-peak elements form an 
almost 56 Ni-empty hole (e.g.. INomoto et al. 1984|> . 

In the sub-Chandra model, the ignition density is too low for electron capture to 
take place. The neutron excess is produced only by the initial CNO elements which are 
converted to 14 N and to 22 Ne, thus depending on the initial metallicity. As a result, the 
mass fraction of 58 Ni is as small as ~ 0.01 (e.g., |Shigcya ma~ et al. 1992] ). 

Such a difference in the amount of 58 Ni can be observationally investigated by late- 
phase (~ lyr since the explosion) spectroscopy at near-infrared (NIR) wavelength. Be- 
cause the ejecta become optically thin in late phases, spectroscopy provides an unbiased, 
direct view of the innermost regions. 

Figure [3] shows the spectral feature around ~ 7000 — 7500 A, i.e., [Fe II] A7155 (with 
some contribution from [Fc II] A7171) and [Ni II] A7378, for 12 SNe la (jMaeda et al. 2010aP 
The [Ni II] A7378 line is emitted from the electron capture region of the ejecta, which is 
supported by the relatively narrow width (< 3, 000 km s _1 ) of the [Ni II] line. Thus the 
existence of [Ni II] line implies the ignition at high density, thus favoring the Chandra 
model. 

It is also interesting to note that these emission lines show the velocity shift, which indi- 
cates the off-center ignition and aspherical nature of SN la explosions ( Maed a~et al. 2010a[ 
2010b). The aspherical features could be seen in the light echo from Tycho's supernova 
remnant (jUsuda et al. 201 i| . 



7. Carbon in SNe la 



Parrent et al. (2011) and Folatelli et al. (2012) investigated the presence of unburned 



material in early-time spectra of SNe la. They find that at least 30% of the objects in 
the sample show absorption at about 6300 A which can be associated with C II A6580 
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(Figure 0]) • This would imply a larger incidence of carbon in SN la ejecta than previously 
noted. If confirmed as carbon, the material producing the observed features would be 
present at very low expansion velocities, merely ~1,000 km s _1 above the Si II velocities. 

Carbon must be present in very deep regions, corresponding to velocities as low asi>« 
11,000 km s _1 . This is well below the expected limit imposed by one-dimensional models, 
and points directly to large mixing effects and/or possible departures from spherical 
symmetry or dumpiness. 

In spherically symmetric models, irrespective of the details of the flame propaga- 
tion (deflagration or detonation), the production of a large amount of 56 Ni (~ 0.6M Q ) 
requires that the strength of the flame, which must lead to the total consumption 
of carbon below 15,000 - 20,000 km s" 1 (jNomoto et al. 19841 |Shigeyama et al. 1992} 
llwamoto et al. 1999[) . The situation is different for non-spherical models. For example, 
the off-center ignition model (IMaeda et al. 2010a| see also lKasen et al. 2009[) synthesizes 
0.54M Q of 56 Ni and carbon with a mass fraction of ~ 0.1 at velocities as low as ^13,000 
km s -1 . This is still larger than the observed velocity, but suggests that the non-spherical 
effects may be important to understand the detection of carbon deep in the ejecta. 

The existence of carbon at relatively low velocity suggests the presence of fair amount 
of unburned C+O materials. The amount Fe peak elements in such an unburned layer de- 
pends on the metallicity. Further observations of early UV spectra could show significant 
metallicity effects. 



8. Concluding Remarks 

As summarized above, recent modeling shows that DD merging could result in both 
the Chandra and the sub-Chandra explosions depending on whether a carbon detonation 
is induced near the surface of more massive WD. The SD scenario could also result in 
both the Chandra and sub-Chandra explosions depending on whether the He shell-flashes 
near the surface of the WD induce a He detonation. 

• For DD-subCh, it is critically important to confirm the formation of surface carbon 
detonation by means of 3D hydrodynamical simulations rather than SPH method. 

• For DD-Ch, whether carbon ignition can occur for both Case 1 and Case 2 needs 
further study for various set of (Mi,M 2 ). 

• For SD-subCh, a mechanism to avoid the production of too much 56 Ni in the surface 
should be studied. In case of the He- WD and C-WD merger, formation of an extended 
He envelope needs to be avoided. 

• For SD-Ch, the outcome of quiescent He-shell flashes, e.g., nucleosynthesis, the rate 
of the He wind mass loss, needs to be studied. 

• Finally, detailed hydrodynamical modeling of the WD spin-down is necessary. In 
this spin-down scenario (e.g., Uustham 20111 IDi Stefano et al. 20111 lllkov fe Soker 20121 
Hac hisu et al. 20T2|) . an almost uniformly rotating C+O WD with a mass range of 1.38 
- 1.48 Mq forms and eventually contract to ignite carbon after a long spin-down time. 
The outcome depends on the ignition density. If it is as high as ~ 10 10 g cm~ 3 , electron 
capture induces collapse rather than explosion (jNomoto fe Kondo 199 If) . If it is lower, 
an SN la explosion would result. 

This research has been supported in part by World Premier International Research 
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JSPS (23540262) and MEXT (22012003, 23105705). 
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